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ABSTRACT 

Context. Water is a key constituent of star-forming matter, but the origin of its line emission and absorption during high-mass star formation is not 
well understood. 

Aims. We study the velocity profiles of low-excitation H2O lines toward 19 high-mass star-forming regions and search for trends with luminosity, 
mass, and evolutionary stage. 

Methods. We decompose high-resolution Herschel-iSFI line spectra near 990, 1110 and 1670 GHz into three distinct physical components. Dense 
cores (protostellar envelopes) are usually seen as narrow absorptions in the H 2 1113 and 1669 GHz ground-state lines, the H 2 987 GHz excited- 
state line, and the H^O 1102 GHz ground-state line. In a few sources, the envelopes appear in emission in some or all studied lines, indicating 
higher temperatures or densities. Broader features due to outflows are usually seen in absorption in the H 2 1113 and 1669 GHz lines, in 987 GHz 
emission, and not seen in Hi 8 0, indicating a lower column density and a higher excitation temperature than the envelope component. A few 
outflows are detected in H' 8 0, indicating higher column densities of shocked gas. In addition, the H 2 1113 and 1669 GHz spectra show narrow 
absorptions by foreground clouds along the line of sight. The lack of corresponding features in the 987 GHz and H^O lines indicates a low column 
density and a low excitation temperature for these clouds, although their derived H 2 ortho/para ratios are close to 3. 

Results. The intensity of the ground state lines of H 2 at 1113 and 1669 GHz does not show significant trends with source luminosity, envelope 
mass, or evolutionary state. In contrast, the flux in the excited-state 987 GHz line appears correlated with luminosity and the H, s O line flux appears 
correlated with the envelope mass. Furthermore, appearance of the envelope in absorption in the 987 GHz and H' 2 g O lines seems to be a sign of an 
early evolutionary stage, as probed by the mid-infrared brightness and the L bot /M em ratio of the source. 

Conclusions. The ground state transitions of H 2 trace the outer parts of the envelopes, so that the effects of star formation are mostly noticeable 
in the outflow wings. These lines are heavily affected by absorption, so that line ratios of H 2 involving the ground states must be treated with 
caution, especially if multiple clouds are superposed as in the extragalactic case. The isotopic H^O line appears to trace the mass of the protostellar 
envelope, indicating that the average H 2 abundance in high-mass protostellar envelopes does not change much with time. The excited state line 
at 987 GHz increases in flux with luminosity and appears to be a good tracer of the mean weighted dust temperature of the source, which may 
explain why it is readily seen in distant galaxies. 

Key words. ISM: molecules - Stars: formation - astrochemistry 



1. Introduction 

The water molecule is a key constituent of star-forming mat- 
ter with a great influence on the formation of stars and plan- 
eto In the gas phase, water acts as a coolant of collapsing 
dense interstellar clouds; in the solid state, it enhances the co- 
agulation of dust grains in protoplanetary disks to make plan- 
etesimals; and as a liquid, it acts as a solvent bringing organic 
molecules together on planetary surfaces, which is a key step 
towards biogenic activity. The first role is especially important 
for high-mass star formation which depends on the balance be- 
tween the collapse of a massive gas cloud and its fragmentation 
(Zinnecker & Yorke 2007). This balance depends strongly on the 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA 

1 This paper uses "water" to denote the chemical species, and H 2 0, 
H^O and HDO to denote specific isotopologues. Unless otherwise in- 
dicated, the gas phase is meant. 



temperature (through the Jeans mass), and the sensitivity of the 
H2O abundance to the temperature, much larger than for CO, 
should make it a useful probe of the high-mass star formation 
process, which is mostly unexplored due to observational diffi- 
culties. 

Interstellar H2O is well known from ground-b ased observa- 
tions of maser emission at c entimeter (22 GH z; Cheung et al.l 
119691) . millimeter (183 GHz: ICernicharo et al.lll990l) . and sub- 
millimeter (325 GHz; Ment en etal.1 Il990h wavelengths. The 
high intrinsic brightness of maser emission makes it useful as 
a signpost of dense g as and for kinematic studies of protostel - 
lar environments (e.g., [Trinidad et aLl l2003 ; San na et al.l|20 12a). 
With VLBI techniques, the proper motions and parallaxes of 
H2O masers can be measured to micro-arcsecond accuracy, lead- 
ing to accurate dis tance estima tes for star-forming regions as 
distant as ~10kpc (San na et alJ 2012b) and a revised picture of 
Galactic structure (Reid et al. 2009)). 



Thermal H2O lines are useful as probes of physical condi- 
tions and the chemical evolution of star-forming regions, but 
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Table 1. Source sample. 



Source 



RA (J2000.0) 
h m s 



Dec 

a I II 



Vlsr 
kms~' 






d 
kpc 



Ref. 



mld-IR-quiet HMPOs 


















IRAS 05358+3543 


05 39 13.1 


+35 45 50 


-17.6 


6.3 


X 


10 3 


1.8 


1. 1 


IRAS 16272^1837 


16 30 58.7 


-48 43 55 


-46.2 


2.4 


X 


10 4 


3.4 


2,2 


NGC 6334I(N) 


17 20 55.2 


-35 45 04 


-4.5 


1.9 


X 


10 3 


1.7 


3,4 


W43 MM1 


18 47 47.0 


-01 54 28 


+98.8 


2.3 


X 


10 4 


5.5 


5,6 


DR21(OH) 


20 39 00.8 


+42 22 48 


-4.5 


1.3 


X 


10 4 


1.5 


7,8 


mid-IR-brlght HMPOs 


















W3 IRS5 


02 25 40.6 


+62 05 51 


-38.4 


1.7 


X 


10 5 


2.0 


9, 10 


IRAS 18089-1732 


18 11 51.5 


-17 3129 


+33.8 


1.3 


X 


10 4 


2.3 


11, 12 


W33A 


18 14 39.1 


-17 52 07 


+37.5 


4.4 


X 


10 4 


2.4 


11,30 


IRAS 18151-1208 


18 17 58.0 


-12 07 27 


+32.0 


2.0 


X 


10 4 


2.9 


13, 13 


AFGL 2591 


20 29 24.7 


+40 11 19 


-5.5 


2.2 


X 


10 5 


3.3 


14,8 


Hot Molecular Cores 


















G327-0.6 


15 53 08.8 


-54 37 01 


-45.0 


5.0 


X 


10 4 


3.3 


15, 16 


NGC 63341 


17 20 53.3 


-35 47 00 


-7.7 


2.6 


X 


10 5 


1.7 


3,4 


G29.96-0.02 


18 46 03.8 


-02 39 22 


+98.7 


3.5 


X 


10 5 


6.0 


17, 18 


G31.41+0.31 


18 47 34.3 


-01 12 46 


+98.8 


2.3 


X 


10 5 


7.9 


19,20 


Ultracompact HII Regions 


















G5.89-0.39 (W28A) 


18 00 30.4 


-24 04 02 


+ 10.0 


5.1 


X 


10 4 


1.3 


21,22 


G10.47+0.03 


18 08 38.2 


-19 5150 


+67.0 


3.7 


X 


10 5 


5.8 


23,20 


G34.26+0.15 


18 53 18.6 


+01 14 58 


+57.2 


3.2 


X 


10 5 


3.3 


24,25 


W51N-el 


19 23 43.8 


+ 14 30 26 


+59.5 


1.0 


X 


10 5 


5.1 


26,27 


NGC7538-IRS1 


23 13 45.3 


+61 28 10 


-57.4 


1.3 


X 


10 5 


2.7 


28,29 



Notes. The first reference is for the luminosity, the second for the distance. If the distance measurement is more recent than the luminosity estimate, 

the luminosity has been scaled to the new distance. The text uses "short" source names, which is the part preceding th e + or - sign. 

References: (1) Beutheretal. (2002)jJ2) foamy et al.1 ((2 002); (3) Sandell (2000); (4) Neckel ( 1 9781); (5) [Nguyen Luong et all feOllh: (6 ) 



lLesteretal.llll985 '): (7) Ja kob et al.l (120 07): (8) Rvgl et al.ld2oTa): (9)lLadd et al.Ul993l) : (10) Hachis uka et alJ <2006h ; (1 DlFaundez etj 



(12)IXuetal.l < l2011l) ; (13) Sridhara rTeTdl (120021) ; (14) Van der Ta ket all ( 119991): (15) [U rquhart et al. (2012); (16) iMinier et all (T2009); (17) 
Cesaroni et al] £l998): (1 8) iPratapet al.l (11 999): (19) M ueller et al.| ( 2002): (20) Churchwell et al. ( 1990W21) IVan der Tak etal] 72 000) (22) 
Mot ogietaljj2011l) ;(23) Hunt er et all T200 0); (24) Hatchell & van der Tak (2003) (25) Kuchar & Bania ( 1994); (26) IVan D ishoeck e t aljd2011l) 
(27) IXu et all J20091) : (28)lSandell & Sievers (2004); (29) Moscadelli et al. (2009); (30) Im mer et al.ld2013l) . 



generally cannot be observed from the ground. Before Herschel, 
space-based submm and far-IR observat ions of H ?,0 lines 
were made with ISO ([V an Dishoeck & Helmich 1996), SWAS 
dMelnick & Ber gin 2005} and Odin (JBjerkeli et al. 2009), but 
these data do not have sufficient angular resolution to determine 
the spatial distribution of H2O. In contrast, space-based mid- 
IR and ground-based mm-wave observations of thermal H2O 
and H^O lines have high angular resolution but only probe the 
small f raction of the gas at h i gh tem peratures (Van der Tak et al. 
2006t IWatson et all 120071; Lf0rgensen &~va n Dishoeck 2010; 
Wangetal.1120121) . 

The first results from the Herschel mission have demon- 
strated the potential of H2O observations of high-mass star- 
forming regions at high spatial and spectral resolution. Mapping 
of the DR21 region shows orders of magnitude variations in H2O 
abundance between various physical components (envelope, out- 
flow, foreground clouds) due to freeze-out, evaporation, warm 
gas-phase chemistry, and photodissociation dVan der Tak et al.l 
2010). Multi-line observations of the high-mass protostar W3 
IRS5 show broad emission as well as blueshifted absorption 
by H2O in the outflow, in particular in lines from excited 
states of H2O, jinderlining the importance of shock chem- 
istry for H2O (Chavarriaet al. 2010). Observations of the p- 
H2O and p-H^O ground state lines and the P-H2O first ex- 
cited state line at 987 GHz toward four high-mass star-forming 



regions indicate H2O abundances ranging from 5x10 to 
4xl0~ 8 witho ut a clear trend with the physical properties of 
the sources ([Marseille et al.l 1201 Obi) . Finally, observations of 
multiple H2O, Hj 8 and Hi O lines toward the massive star- 
forming region NGC 63341 show high-excitation line emission 
from hot (~200K) gas, and a H2O abundance ranging from 
»10~ 8 in cold quiescent g as to 4xl0~ 5 in warm outflow material 
(Emnrech tinger et aDl2010l) . However, these papers study either 
single sources or small sets of sources, which makes a trend anal- 
ysis inconclusive, if not impossible. 

This paper presents observations of low-excitation lines of 
H2O and H^O toward 19 regions of high-mass star formation 
(Table[TJl. The velocity-resolved spectra show a mixture of emis- 
sion and absorption due to various physical components, which 
we disentangle by fitting Gaussian profiles. The resulting line 
fluxes for the protostellar envelopes and the molecular outflows 
are compared with basic source parameters such as luminosity 
and mass, as well as to the evolutionary stage of the source as 
estimated in various ways. 

To probe the bulk of the material in the protostellar environ- 
ment, this paper focuses on the lowest excitation lines of H2O. 
The advantage of P-H2O is that its two lowest-excitation lines 
are close in frequency (Table |2]l so that the telescope beam is 
similar in size. This coincidence allows us to make estimates 
of the column density and the excitation temperature of H2O 
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Table 2. Observed lines. 



Table 3. SSB continuum flux densities (Jy/beam). 



Molecule Transition Frequency 
GHz 



F 
K 



$mb 0"mb 

arcsec mK 



H 2 


1 1 1 — Ooo 


1113.343 


53.4 


19.2 


40 


H> 8 


1 1 1 — Ooo 


1101.698 


53.4 


19.2 


40 


H 2 


2o2-lii 


987.927 


100.8 


21.5 


50 


H 2 


2l2~ lot 


1669.905 


114.4 


12.7 


75 



Notes. The rms noise level in the last column is for HRS velocity reso- 
lution (0.5 MHz). 



which are almost independent of assumptions on the size or the 
shape of the object. In addition, we use observations of the o- 
H2O ground-state line at 1669 GHz to probe the kinematics of 
the sources, which have a high continuum brightness at this high 
frequency, so that absorption lines are easily detected. This paper 
does not discuss observations of the other 0-H2O ground-state 
line at 557 GHz, because of the large difference in frequency and 
beam size. Maps of the 557 GHz line toward our sources will be 
presented elsewhere. 

Table Q] presents the sour ces, which have been selected 
from several survey s (Molinari et al. 1996 Sridharan et al.l2002l 
IWood & Churchwellll989tlVan der Tak et al . 2000), where pref- 
erence is given to nearby (;2kpc) objects and 'clean' ob- 
jects where one source dominates the emission on «30" scales, 
which is the relevant scale for our observations. Furthermore, 
the sources were selected to cover a range of evolutionary 
stages. In high-mass protostellar objects (HMPOs), the cen- 
tral star is surrounded by a massive envelope with a centrally 
peaked temperature and density distribution. These objects show 
signs of active star formation such as outflows and masers. 
This paper distinguishes mid-IR-bright and -quiet HMPOs, with 
a boundary of 10 Jy at 12 \xva, which may be an evolu- 
tionary difference dVan der Tak et all l2000t iMotte et al.l l2007t 
Lopez-Sepulc re et al.ll2010l) . The presence of a 'hot core' (sub- 
millimeter emission lines from complex organic molecules) or 
an 'ultracompact HII region' (free-free continuum emission of 
more than a few mJy at a distance of a few kpc) are thought to be 
signposts of advanced stages of protostellar evolution. Hot cores 
occur when a star heats its surroundings such that ice mantles 
evaporate off dust grains which alters the chemical composition 
of the circumstellar gas. When the stellar atmosphere becomes 
hot enough for the production of significant ultraviolet radiation, 
the circumstellar gas is ionized and an ultracompact HII region 
is created. 

This paper is organized as follows: Section 2 describes the 
observations and Section 3 their results. Section 4 presents a 
search for trends in the results with luminosity, mass or age, 
and Section 5 provides the discussion and our conclusions. In 
addition to this overview paper, our group is preparing detailed 
multi-line studies of individual sources. 



2. Observations and data reduction 

The sources were obser ved with the Heterodyn e Instrument for 
the Far-Infrared (HIFI; iDe Graauw et al.ll2010l) o nboard ESA's 
Herschel Space Observatory (Pilbratt et al. 2010) in the course 
of 2010 and 2011. Spectra were taken in double sideband mode 
using receiver bands 4b, 4a and 6b, respectively, for the 1 n-Ooo, 
2o2— 1 11 and 2i2— lot lines (Table |2). The lu-Ooo lines of H2O 
and H^O were observed simultaneously. All data were taken us- 
ing the double beam switch observing mode with a throw of 2. '5 



Source 



987 GHz 
(304 yum) 



1108 GHz 
(271 yum) 



1669 GHz 
(180 pm) 



IRAS 05358 


146 


235 


447 


IRAS 16272 


382 


579 


787 


NGC 6334 I(N) 


892 


1100 


887 


W43 MM1 


772 


995 


1073 


DR21(OH) 


896 


1248 


1678 


W3 IRS5 


546 


760 


1896 


IRAS 18089 


393 


507 


837 


W33A 


415 


593 


1216 


IRAS 18151 


106 


217 


308 


AFGL2591 


310 


467 


830 


G327-0.6 


1206 


1426 


1621 


NGC 63341 


2186 


2601 


4624 


G29.96 


510 


688 


1372 


G31.41 


492 


684 


948 


G5.89 


1020 


1359 


2908 


G10.47 


914 


1340 


2368 


G34.26 


1566 


2272 


3549 


W51N-el 


1894 


2388 


2468 


NGC 7538 IRS 1 


474 


633 


1037 



Notes. The values for 1 108 GHz are averages of measurements near the 
H 2 1113 GHz and the H; 8 lines. 



to the S W, as part of the Guaranteed Time k ey program Water in 
Star-f orming regions with Herschel (WISH: [Van Dishoeck et aD 
|2011|) . See Appendix lAl for a detailed observing log. 

Data were simultaneously taken with the acousto-optical 
Wide-Band Spectrometer (WBS) and the correlator-based High- 
Resolution Spectrometer (HRS), in both horizontal and vertical 
polarization. This paper focuses on the HRS data, which cover 
230 MHz bandwidth at 0.48 MHz (-O.Bkms 1 ) resolution, al- 
though for a few sources where the 1669 GHz line is too broad 
to fit in the HRS backend, the WBS data are used, which cover 
1140MHz bandwidth at 1.1 MHz (~0.30kms _1 ) resolution. 
Spectroscopic parameters of the lines are taken from the JPL 
(Pickett etal. 1998) and CDMS (Muller e t al. 2001 ) catalogs 
The adopted FWHM beam sizes are taken from Roelf sema et al.l 
(12012,) and scaled to our observing frequencies. System tem- 
peratures are 340-360 K (DSB) at 1113 GHz, 400^150 K at 
988GHz, and 1400-1500K at 1670GHz. Integration times are 
6.4 minutes at 1113 GHz (ON+OFF), 12 minutes at 988 GHz 
and 34 minutes at 1670 GHz. 

The calibration of the data was performed in the Herschel 
Interactive Processing Environment (HIPE) version 6 or higher; 
further analysis was done within the CLASSo package, version 
of January 201 1. The intensity sca le was converted t o T m t, scale 
using main beam efficiencies from Roelfse ma et al.1 (1201 2l) . and 
linear baselines were subtracted. After inspection, data from the 
two polarizations were averaged together to obtain rms noise 
levels reported in the last column of Table [2] The calibration 
uncertainty is 15% at -1000 GHz and 20% at -1700 GHz (V. 
Ossenkopf and R. Shipman, priv. comm.fl 
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Spectral Index 

Fig. 1. Spectral indices of our sources, calculated from the 
HIFI continuum levels at 987 and 1669 GHz, versus the mass- 
weighted temperatures calculated from the models in § 14.11 



3. Results 

3.1. Continuum emission 

We have measured the continuum levels in the HIFI data by fit- 
ting zeroth- or first-order polynomials to the velocity ranges in 
the spectra without any appreciable line signal. In the cases of 
G327-0.6, NGC 63341, and NGC6334I(N), the HRS backend is 
not broad enough to contain the 1669 GHz line, so that we have 
used the WBS data instead. In general, the differences between 
continuum levels from the two backends (and within polarization 
channels from a given backend) are within 5% near 1000 GHz 
and within 10% at 1670 GHz, which is within the calibration 
uncertainty (©. The flux densities S v in Jy/beam were derived 
from the measured DSB continuum levels by 

„ * kT A 



2nr, a D 2 



dRohlfs & Wilsonll2004 Eq. 7.18) where the effective telescope 
diam eter D - 3 . 28 m and the aperture efficiency r] a =s 0.65 
(Roelfsemaet al. 2012) with a small frequency dependence 
due to Ruze losses. The factor 2 converts DSB to SSB sig- 
nal, and T A — TjfT^ with the forward efficiency r\f - 0.96 
dRoelfsema et al.ll2012l) . 

As seen from Table [3] the continuum flux densities of our 
sources generally increase with frequency, as expected for ther- 
mal emission from warm dust. The fre e-free emission of our 
sourc es is ! 1 Jy which is negligible dVan der Tak & Mentenl 
2005). For optically thick dust emission at a temperature Tj, 
S v - B v (T,i) so that in the Rayleigh-Jeans limit, the spectral 
index a (defined as S v oc y°) equals 2, while for optically thin 
emission, the spectrum steepens to a — 2 + j3 — 3^1, depend- 
ing on the far-infrared dust opacity index /?, which is 1-2 in 
most dust models (Ossenkopf & Henning 1994). However, at the 
high frequencies of our observations and the low average tem- 
peratures of our sources (15-40 K: Table |4), deviations from 
the Rayleigh-Jeans approximation are significant, which reduce 
these expected spectral indices by «1 at Tj — 40 K and by «2 at 
T d = 15 K. 




-60 -40 

Velocity (km/s) 

Fig. 2. Gaussian decomposition of the observed H2O line pro- 
files (after continuum subtraction) toward IRAS 16272. The en- 
velope component is drawn in red, the broad outflow in green, 
the narrow outflow in blue, and foreground clouds in purple. 



The spectral indices of our sources, calculated between 987 
and 1670 GHz, are mostly between 1.2 and 2.4 (Fig.[TJ, which 
indicates that the emission is close to optically thick. The mod- 
els in 94.1l also indicate continuum optical depths around unity at 
-1000 GHz. Four sources (NGC 6334I(N), W43 MM1, G327- 
0.6, and W51e) have lower spectral indices (between 0.0 and 
0.5), which for compact sources would suggest low dust tem- 
peratures. For extended sources (size ! 15"), the emission would 
be resolved at the highest frequencies, which would artificially 
flatten the spectrum. The model calculations in JH.ll indicate me- 
dian source sizes of 17.6, 14.2, and 8.5" FWHM at 987, 1113, 
and 1669 GHz, which is 83%, 75%, and 67% (i.e., most) of the 
corresponding beam size, suggesting that size effects play a mi- 
nor role. The importance of dust temperature is also suggested 
by Fig. Q] which shows that the sources with low spectral in- 
dices have low mass-averaged temperatures as calculated from 
the models in 94.11 We conclude that our observed variations in 
the spectral index of the HIFI continuum emission are not due to 
variations in the source size, but rather to a low dust temperature, 
both directly and through deviation from the Rayleigh-Jeans law. 



3.2. H 2 lines 



2 [http://www.iram.fr/IRAMFR/GILDAS] 

3 See http://herschel.esac.esa.int/twiki/bin/view/ 
Public/HifiCalibrationWeb?template=viewprint 



Figures IB.U4B.19l in Appendix [B] show the observed line pro- 
files. All four lines are detected toward all nineteen sources (ex- 
cept H^O toward IRAS 18151), and the line profiles are a mix- 
ture of emission and absorption features. Analogous to the case 
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Fig. 3. Central velocities of the envelope, narrow outflow, and 
broad outflow components of the observed H2O lines compared 
with each other and with the systemic velocity from ground- 
based data. In the bottom right-hand comer, the typical error bar 
on the velocities is shown. 



of DR21 dVan der Tak et al.l 12010). three physical components 
appear from these observations (Fig. |2). 

First are the protostellar envelopes (also sometimes called 
dense cores) which usually appear in all four lines in absorption 
at an LSR velocity known from ground-based millimeter-wave 
emission line data (Fig. [3] top right); the widths are also similar 
to the ground-based values. Several sources show two envelope 
components, due to binarity (e.g., W3 IRS5) or infall/outflow 
profiles (e.g., W33A, G34.26). The general appearance in ab- 
sorption indicates a high H2O column density and a low excita- 
tion temperature for this component, which is further discussed 

in £OD 

Second are the molecular outflows which usually appear 
in 11 13 and 1669 GHz absorption and 987 GHz emission at 
velocity offsets of 5-10 km s~' from the envelope signal and 
with widths of 10-20 km s _I FWHM or sometimes more, up 
to 40kms _I . Most sources show two outflow components: a 
broad (AV ~20kms~ I ) and a narrower (AV ~10kms _1 ) fea- 
ture. In some cases, the width of the narrow outflow feature 
is similar to that of the corresponding envelope, but the com- 
parison to the ground-based velocity and the presence or ab- 
sence of Hj 8 emission or absorption allows for an unambigu- 
ous assignment. The lack of H^O signal and the appearance in 
987 GHz emission suggest that the outflows have lower H2O col- 
umn densities and higher excitation temperatures than the en- 
velopes. Section [3~5l discusses the outflow component in detail. 

Third are foreground clouds, which appear as narrow 
(AV<5kms _I ) 11 13 and 1669 GHz absorption features at veloc- 
ities offset by 5-50kms _I from the envelope signal. The lack 
of signal in the 987 GHz and Hi O lines indicates a low H2O 
column density and a low excitation temperature for these com- 
ponents, without exception. The number of foreground clouds 
within the HRS bandwidth ranges from zero toward IRAS 05358 
and NGC 7538 to seven toward W51e. Many of these fore- 
ground clouds are known from ground-based mm-wave ob- 
servations of ground- sta te molecular lines such as CS, HCO + 
and HCN dGreaves & Williamsl [l994; Godard et al. 2010b and 
from [HI] 21 cm observations dFish et all i2003; Pandian et al.l 




AV (envelope H 2 0) [km/s] 

Fig. 4. Line widths (FWHM) of the envelope components of the 
observed H2O lines c ompared to the widths of the C 18 9-8 line 
measured with HIFI dSan Jose Garcia et al.l l2.Q13). Error bars are 
smaller than the plotting symbols. 



2008). Other foreground clouds have first been seen with HIFI, 
such as the V=+24kms~ 1 c loud toward W51e which is seen 
in HF dSonnentrucker et all [2010) but not in our H2O data, 
and the V=+13kms _1 cloud towar d AFGL 2591 which i s seen 
in H2O and HF but not in CO (Emrjrechtingeret al. 2012; 
IVan der Wiel et al.1 120131; IChoiet all 120131). Some foreground 
clouds are intervening objects on the line of sight, while others 
represent the 'outer envelopes' of large-scale molecular cloud 
complexes, such as the V=+65kms _1 cloud toward W51e dKoo 
1997) and two clouds toward NGC 6334 dVan der Wiel et al 
2010). See §13.6l for further discussion of the foreground clouds. 
Our decomposition of the line profiles is similar to the 
description of H2O line profiles toward low-mass protostars 
in terms of 'narrow' , 'medium' and 'broad' components 
dKristensen et alJl2010h . where narrow = envelope, medium = 
narrow outflow, and broad = broad outflow. Because of the large 
range in luminosity and mass within our sample, we base the as- 
signment of features to envelope or outflow not only on the line 
width, but also on the appearance in emission or absorption in 
the various lines. 



3.3. Profile decomposition and flux extraction 

We have measured the fluxes (jT m bdV; in emission) and ab- 
sorbances ( f rdV; in absorption) of the lines by fitting Gaussians 
to the observed profiles for each line independently. Figure [2] 
shows IRAS 16272 as an example of the procedure, and the re- 
sults are presented in Table lA.2l for the envelopes, Table IA.3~l for 
the outflows, and Table lA~4l for the foreground clouds. Our fit- 
ting procedure is similar to that by Flagev et al. (2013), where 
the reader can find further details. 

Since some of the emission and absorption features in our 
data have non-Gaussian shapes, we have verified the results of 
the Gaussian fitting procedure by a "moment" analysis, which 
just adds up the flux in all spectrometer channels within a given 
velocity range. This alternative method extracts the flux of a 
spectral component independent of its shape, but it cannot han- 
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die blending of components, which occurs frequently in our 
data, unless further assumptions are made. Another disadvan- 
tage of the moment method is that the definitions of the veloc- 
ity ranges are often somewhat arbitrary. After some experiment- 
ing, we conclude that fitting Gaussians gives more reliable re- 
sults than taking moments, if we accept that in some cases, one 
physical component (usually the outflow) is represented by sev- 
eral (up to three) Gaussians. Another consequence of the non- 
Gaussian line shape is that the extracted central position and the 
width of a given component vary somewhat from line to line: up 
to 2kms~ 1 for the envelopes and up to 5kms -1 for the broad 



outflow component. The uncertainties on the line positions and 
widths in Tables lA2llA.3l and lA.4l represent this line-to-line vari- 
ation, whereas the uncertainties on the fluxes and absorbances 
are dominated by the overall calibration uncertainty (©, which 
is always larger than the formal error of the Gaussian fit caused 
by the deviation of the observed profile from a Gaussian shape. 
For absorption features, Tables [A2~l4A.4| give the absorbance 
(= the velocity-integrated optical depth): 



/ 



rdV = -1.06AVln( c — ' ) 



where 7) is the peak depth of the feature, T c is the SSB con- 
tinuum temperature, and AV is the FWHM line width. For small 
optical depth, the error on 7/ propagates linearly into the ab- 
sorbance, but if the absorption is almost saturated (7) — > T c ), the 
error on the absorbance is much larger than that on the measured 
absorption depth. Cases where the error exceeds the absorbance 
due to saturation, which occurs mostly for envelope components, 
are marked with lower limits in Tables IA.21 - lA.4l 

In the frequent case that absorption and emission features 
are partially blended, we have estimated the line contribution 
to the absorption background from the spectra, and added the 
result to the continuum brightness. This procedure assumes that 
the absorber is located in front of both the continuum and the 
line emitting source, which is plausible since in most cases, the 
absorption depth significantly exceeds the continuum brightness 
alone. 

The model calculations in JH.ll indicate median source sizes 
of 17.6, 14.2, and 8.5 arc seconds FWHM at 987, 1113, and 
1669 GHz, which is 83%, 75%, and 67% of the corresponding 
beam size. The filling factor would be the square of this fraction 
which is 68%, 56%, and 45% respectively, so that our absorp- 
tion line measurements will be affected by covering factor ef- 
fects if the absorbers are much smaller in size than the beam. In 
this case, the beam-averaged optical depth would be below the 
source-averaged value. 

3.4. Envelopes 

The H2O line spectra of all sources show envelope components, 
and the velocities and widths of these components in Table I A. 21 
agree well with previous measurements for these sources in other 
molecules from the ground (Fig. [3}. The widths also mostly 
agree with the widths of the C 18 9-8 line observed with HIFI 
in a similar beam (Fig. 0J; the few cases where AV(C 18 0) > 
AV(H20) may be influenced by an outflow contribution to the 
C 18 9-8 line. The measured H2O line widths range over a fac- 
tor of 4, and show a clear correlation with the envelope masses 
calculated in j4~T1 (Fig.l5Tl. 

The appearance of the envelopes in the four lines differs 
markedly from source to source, which can be grouped into 
four types. Seven envelopes appear in absorption in all four 
lines: G31.41, DR210H, G327, NGC6334I(N), IRAS 05358, 
W43MM1, and one of the G34.26 envelopes. Six envelopes 
show emission in the 987 GHz line and absorption in the other 
lines: G10.47, IRAS 18089, W51e, NGC 63341, IRAS 16272, 
and one of the W33A envelopes. A further six envelopes appear 
in H 2 987 GHz and in H^ 8 emission, and in H 2 1113 and 
1669 GHz absorption: G29.96, NGC 7538, W3 IRS5, G5.89, 
and the second envelopes of W33A and G34.26. Finally, one 
source, AFGL 2591, appears in emission in all four lines. One 
curious case is IRAS 18151, where the non-detection of the 
Hl 8 line leaves ambiguous whether it belongs in the second 
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Fig. 6. Top: H^O line profile toward G5.89, showing a P Cygni 
profile. Bottom: H^O line profile toward G10.47, showing an 
inverse P Cygni profile. The bottom spectrum is continuum sub- 
tracted. Vertical scale is r m b in K, and the dotted line denotes the 
systemic velocity from ground-based data. 



or third of the above source types; or perhaps it is in a transition 
from the second to the third type. 

One possible origin for the diverse appearance of the enve- 
lope features is related to temperature, as appearance in emis- 
sion (absorption) requires an excitation temperature above (be- 
low) the continuum brightness. In this scenario, pure absorption 
sources would have the most massive and coldest envelopes, 
sources where the 987 GHz line appears in emission would be 
somewhat warmer, sources where both H2O 987 GHz and Hi 8 
appear in emission would be even warmer, and the pure emission 
source would be the warmest. Alternatively, the appearance of 
the envelope features in emission or absorption may be a matter 
of continuum brightness, as detectability in absorption requires a 
background source whereas emission does not. In particular, the 
presence of outflow cavities in the HIFI beam would decrease the 
continuum compared with a beam filled with envelope emission. 
However, the continuum flux densities in Table [3] do not sup- 
port this option: the four envelope types above have very sim- 
ilar median continuum brightness levels. Outflow cavities may 
also affect the H2O line signals, but the net effect may either be 
emission or absorption, depending on the difference between the 
H2O excitation temperature and the dust temperature, and on the 
relative location of dust and water along the line of sight. We 
conclude that heating effects are more likely than geometric ef- 
fects to be responsible for the variation in envelope appearance 
between our sources. 
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Fig. 7. Line widths (FWHM) of the envelope, narrow outflow, 
and broad outflow components of the observed H2O lines com- 
pared with each other. Error bars are smaller than the plotting 
symbols. 



Some sources show two envelope features in some or all 
H2O line profiles, which in certain cases may be due to the pres- 
ence of common-envelo pe binaries at ~ 1000 AU separation, as 
suggested for W3 IRS5 (Chavarriaet al. 2010). In other cases, 
the lines have P Cygni or inverse P Cygni profiles, which are 
due to expansion and infall, respectively (Fig. |6}. Such pro- 
files have also b een observed in the H7 O line spectra of low- 
mass protostars (Kristensen et al. 2012b and for the low-mass 



pre-stellar core L1544 (Casell iet al.l2012l) . Interestingly, the (in- 
verse) P Cygni signature is usually not seen in all lines toward 
a given source. The sources G10.47 and G34.26 show inverse 
P Cygni profiles in the H2O 987 GHz and H^O lines, suggesting 
that infall motion only takes place in their inner envelopes. For 
G34.26, this result is confirmed by SOFIA obse rvations of red - 
shifted NH 3 line absorption at 1810 GHz dWvrowski et al.l2012l) . 
and f or G 10.47 by SMA observations of vibrationally excited 
HCN dRoiffsetaLlboill) . In contrast, G5.89 and DR210H show 
P Cygni profiles in the H2O 1113 and 1669 GHz lines and in 
H^O, but not in the H 2 987 GHz line, suggesting that expand- 
ing motions only take place in their outer envelopes. In W33A, 
where only the Hj 8 line may show a P Cygni profile, the ex- 
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AV (broad H 2 0) [km/s] 

Fig. 8. Line widths (FWHM) of the broad outflow components 
of the observed H2O lines compar ed to the widths of the CO 3 - 
2 lines measured with the JCMT (San Jose Garcia et al.ll2013b . 
Error bars are smaller than the plotting symbols. 



panding motion may take place further inside the envelope. For 
W43 MM1, the low-/ lines do not sh ow infall or outflow mo- 
tions, consistent with the observation bv lHerpin et al.l (120 12) that 
only high-/ lines show infall motions in this source, i.e., that in- 
fall is confined to the inner envelope. Similar qualitative kine- 
matic variations with radius have been observed in HCN lines 
toward the Sgr B2 envelope, where the emission peak shifts from 
blue to red with increasing /, suggesting infall in the outer and 
expansion in the inner envelope dRolffs et alJl2010h . 

3.5. Outflows 

The outflows usually show up in the H2O line profiles as 
two components: a broad (FWHM -10-20 km s _1 ) emission 
component visible in all lines, and a narrower (FWHM ~5- 
lOkms 1 ) absorption component seen in all but the 987 GHz 
line. Exceptions are AFGL 2591, IRAS 18151 and IRAS 05358, 
which only show a broad component, and G5.89, which has 
an exceptionally broad and complex line profile; t his source is 
know n to have a very massive and powerful outflow (lAcord et al.l 
fT997h . For the sources DR210H, W3 IRS5 and NGC 7538, the 
widths of the narrower components are so small (2-3kms _1 ) 
that these may be 'second envelopes', as also seen toward W33A 
and G34.26, and as also advocated by IChavarria et al.l (|2010) 
for W3 IRS5. The mm-wave continuum emissio n from these 
sourc e s indeed shows multiple spatial components (Woody et alJ 
1989; IVan der Tak et al.1 20051). sometimes with mi d-infrared 
counterparts (Megeath et al.l l2005b Ibe Wit et al.ll2009h . 

For most of our sources, the broad outflow component ap- 
pears in emission in all four lines, but there are some exceptions. 
For AFGL 2591, NGC 63341, W51e and G34.26, the broad out- 
flow component appears in absorption in the 1113 and 1669 
GHz lines, and in emission in the 987 GHz line; for G 10.47, 
this component appears in emission in the 1113 and 987 GHz 
lines, and in absorption in the 1669 GHz line. Toward W51e, 
a 'very broad' component appears in emission in all four lines 
between +30 and +90 kms -1 . It is unclear in which category 



G31.41 belongs, since the broad outflow component is only de- 
tected in the 987 GHz line, perhaps because in the ground state 
lines, emission and absorption cancel each other out. The shape 
of the broad outflow component on the 987 GHz line profile is 
often non-Gaussian, but the outflow contribution can be modeled 
with two Gaussians of about the same height and width. In this 
case the velocity in Table lA3l is the mean of the values for the 
two Gaussians, while the width and the flux are the sums. 

The narrow outflow component usually appears in absorp- 
tion in the H 2 1 1 13 and 1669 GHz lines and the H' s O line, and 
is not seen in the H2O 987 GHz line, but again, exceptions occur. 
For W3 IRS5, G327 and G31.41, the narrow outflow component 
appears in 1113 and 1669 GHz absorption and in 987 GHz emis- 
sion; for IRAS 18089 and G34.26, this component appears in 
1113 and 1669 GHz absorption and in H^O emission, while for 
G10.47, it appears in emission in all four lines. Unlike the en- 
velopes, the outflow absorptions are usually not saturated, with 
G327 as only exception. 

We have used the measured velocities of the outflow compo- 
nents to test if our sources have a preferred orientation in the sky. 
In particular, for a face-on geometry, the outflow would be di- 
rected toward us, and thus appear blueshifted relative to the enve- 
lope component if seen in absorption. The narrow outflow com- 
ponent is the one that is usually seen in absorption, and Figure [3] 
shows that this component is indeed usually blueshifted relative 
to the envelope. The exceptions are IRAS 16272, DR21(OH), 
G29.96, and W51e, and Fig. [2] illustrates the case of IRAS 
16272, especially for the 1669 GHz line. The widths of these 
redshifted absorption features are too large to be due to infall 
motions, and we suggest that these sources have binaries or other 
deviations from centrosymmetric geometry that allow their re- 
ceding outflow lobes to show in absorption. In addition, Figure|7] 
shows that there are no clear trends between the widths of the en- 
velope and outflow components. 

Figure [8] compares the widths of the broad outflow com- 
ponents to the widths observed in CO 3-2 with the JCMT 
(San Jose Garcia et al. 2013). While sources with broader H2O 
outflows also seem to have broader CO outflows, the width for 
any given source is larger in CO than in H2O by about 10 km s . 
This result is contrary to the case of low -mass protostars, wher e 
the H2O line widths exceed those of CO (Kristen sen et al.l2 012). 
For our sources, the H2O may be confined to denser gas where 
velocities are lower. Alternatively, the difference is excitation: 
our H2O values are averages over ground-state and excited-state 
lines, while the outflow line widths increase with / or E up . 
Indeed, the widths of the excited-state H 2O lines are very close 
to t hose measured in CO , as also found bv lChavarria et al.l (120 1 0) 
and lHerpin et"aTl (120121) . 



3.6. Foreground clouds 

The foreground clouds appear as narrow (FWHM = \^\ km s ) 
absorption features in the 1113 and 1 669 GHz line spectra of our 
sources. While one cloud also shows H^O absorption, the gen- 
eral lack of signal in the H2O 987 GHz and H^ s O lines indicates 
a low column density and a low volume density for these clouds. 
The last two columns of Table|A4]list estimates of the H2O col- 
umn densities and ortho-para ratios of the foreground clouds. 
These values were obtained from the velocity-integrated opti- 
cal depths in the table by assuming that in these clouds, all H2O 
molecules are in their ortho and par a ground states, i.e., that exci - 
tation is negligible. As discussed by Emprechting er et al.l(l2013l) . 
this assumption generally holds for column density estimates 
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Table 4. Parameters of continuum models. 



Input parameters Output parameters 



Source 



"out 

( 10 4 AU) 



V 



(Mo) 



(K) 



IRAS 05358 


5.4 


1.5 


142 


22 


IRAS 16272 


17.0 


1.5 


2170 


18 


NGC 6334I(N) 


7.2 


1.3 


3826 


15 


W43 MM1 


14.7 


1.5 


7550 


21 


DR21(OH) 


5.0 


1.35 


472 


27 


W3 IRS5 


11.3 


1.4 


424 


34 


IRAS 18089 


3.9 


1.9 


172 


41 


W33A 


7.1 


1.2 


700 


26 


IRAS 18151 


5.0 


1.4 


153 


27 


AFGL2591 


7.1 


1.0 


320 


38 


G327 


8.0 


1.2 


2044 


24 


NGC 63341 


6.8 


1.5 


750 


38 


G29.96 


9.6 


1.4 


768 


40 


G31.41 


11.9 


1.4 


2968 


32 


G5.89 


3.6 


1.3 


140 


36 


G10.47 


6.0 


1.4 


1168 


41 


G34.26 


8.1 


1.4 


1792 


39 


W51N-el 


14.0 


1.4 


4530 


38 


NGC 7538 


6.4 


1.0 


433 


36 



Mass- weighted temperature of the envelope model. 



1 E 




10 
Offset [arcsec] 

Fig. 9. Submillimeter brightness profiles of AFGL 2591 mea- 
sured with JCMT/SCUBA (error bars), and synthetic brightness 
profiles derived from WR models and convolved by the SCUBA 
beams for p=l.O (solid line) and p- 1.5 (dashed line). The ver- 
tical dotted line indicates the adopted source radius. For this 
source, we choose a density profile of p=1.0. 



from the 1113 and 1669 GHz lines, but fails for the H 2 557 
GHz line. 

The H2O column densities of the foreground clouds are 
found to lie mostly in the range lxlO 12 - lxlO 13 cm" 2 , with 
a few outliers at the upper an d lower end. These values are com- 
parable to earlier estimates (IVan der Tak et al.l 1201 Ol Lis et alJ 
2010) and indicate that the H2O abundance in these clouds 
is limited by photodissociation. The derived ortho/para ratios 
for H2O lie mostly in the range between 3 and 6, again with 



a few outliers to higher and lower valu es. Our derived o/p 
ratio s are very sim ilar to the values by Lis et al.l (1201 Ol) and 
Herpin et al. (2012). Comparison to the detailed multi-line study 
by Flagev etal] d2013h for the 5 overlapping sources suggests 
that our estimated column densities and o/p ratios for individual 
clouds have uncertainties of a factor of ~2, which is an upper 
limit since not exactly the same positions were observed. 

The ortho-para ratio of H2O is expected to rise from a 1.0 at 
low temp eratures («T5 K) to «3 at high temperatures ( ; 50 K) as 
shown by Mumma et al. ( 1987). Even though the values derived 
here are probably uncertain to a factor of ~2, our results are in- 
consistent with the low-temperature limit. If the H2O molecules 
were formed on cold (10 K) grain surfaces, their spin tempera- 
ture does not reflect these conditions, and may instead be set by 
the grain temperature shortly before the H2O ice mantle evapo- 
rated (at T *100K). Alternatively, the H2O molecules may have 
formed in the warm (7\i n ~50K) gas phase by ion-molecule re- 
actions, or the o/p ratio of the evaporated H2O may have equi- 
librated with that of H? in the warm gas phase in reactions with 
H 3 +. 

4. Analysis and discussion 

4.1. Mass estimation 

Estimates of the envelope masses exist in the literature for many 
of our sources. However, these mass estimations were performed 
by different authors, using different techniques and based on data 
from different telescopes. In order to have a homogeneous mass 
determination over our sample, we derive the envelope mass by 
comparing models of the continuum spectral energy distribution 
(SED) with the available observations for each source. 

We use a modified version of the 3D continuum radiative 
transf er code HOCHUNK3D (IWhifnev et alJ I20H iRobrtaiiiel 
|201 lb (hereafter WR) to find the envelope temperature profile 
and derive a model SED for each s ource. Our modified ve r sion o f 
the code uses dust opacities from Ossenkop f & Henningl dl994). 
model 5 with thin ice mantles an d a gas density of 10 6 c m" 3 , 
suitable for high-mass protostars (lVan der Tak et al.lll999l) . and 
allows to vary the radial density power-law exponent p. The WR 
code allows for a detailed description of the protostellar envi- 
ronment in terms of a star, a disk, a cavity and an envelope. To 
decrease the number of free parameters, we use a simple spher- 
ically symmetric model with no cavity and no disk. This sim- 
plification has its main effect in the near to mid-IR region of the 
spectrum, while the sub-millimeter and millimeter- wave regions, 
where the dust is mainly emitting, will stay almost unchanged. 

The envelope size and density exponent p are derived from 
archival sub-millim eter images from the JCMT/SCUBAo and 
APEX/LABOCA dSiringo et al.ll2009l) instruments. For the en- 
velope size, we choose the radius of the contour at the 3<x noise 
level, while for the density exponent, we use the WR code to 
create synthetic sub-millimeter maps for various p-values. We 
determine the best-fit density exponent by minimizing the differ- 
ence between the observed and synthetic brightness profiles (see 
Fig. |5). The synthetic maps are created using the correspond- 
ing instrumental pass-band and convolved with the instrumental 
beam profile including the error-beam. 

With the outer radius and the density profile fixed, the only 
remaining parameters are the envelope mass and the bolometric 



4 The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 
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Fig. 10. Top: Model (continuous line) and observed SED fluxes 
(error bars) for IRAS 16272. The SED points are normalized by 
the source size for A > 100 microns, as shown by the "X" sym- 
bols. Continuum flux densities from our HIFI data are labeled in 
red. Bottom: temperature and density profiles for IRAS 16272. 



luminosity. The observed fluxes from the literature correspond 
to the peak emission at the source position. In most cases, the 
beam size of the observed fluxes is smaller than the source size. 
To correct for the missing flux, we normalize the observed fluxes 
by the ratio between the flux from the whole model and the flux 
from a region as big as the beam size. As an example, Figure [TOl 
shows the observed and derived SED for the source IRAS 16272 
as well as the temperature profile. Table |4] lists the adopted and 
the derived physical parameters for each source, and AppendixICl 
presents the results in graphical form. 
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and filled purple circles are ultracompact HII regions. 
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Fig. 12. Envelope H2O 987 GHz line flux versus mass-weighted 
temperature in our envelope models. 



4.2. Correlations with luminosity and mass 

To understand the origin of the observed H2O lines, we now in- 
vestigate if the measured line fluxes show trends with the lu- 
minosity of the source Lboi (Table [TJ and the mass of the en- 
velope M env (Table 0). In addition, we have searched for trends 
with distance but did not find any. For the luminosity, we use 
the results of the model fit in 94.11 so that the source flux at all 
wavelengths refers to the same area on the sky. Figures ID.H - 
ID.3I in Appendix [D] show scatter plots of the line fluxes of the 
envelope and the narrow and broad outflow components versus 
Lboi an d M em . In these figures, absorption is plotted as negative 
emission, and 'double envelope' features are plotted as two data 
points. We identify possible trends by visual inspection of these 
plots, and test these (if applicable) by computing the correlation 
coefficient r and the probability of false correlation P, which de- 
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pends on r and the n umber of data points N (see Appendix of 
iMarseille et aTll2010al) . 

For the envelope component, the visual inspection reveals 
two trends. First, the 987 GHz line flux increases with Lboi, but 
only if the line appears in emission. This trend is not a corre- 
lation, but more a 'threshold effect': strong 987 GHz line emis- 
sion is only seen at high luminosities (Lboi - 10 5 L ). Second, the 
H^O line flux decreases with increasing M env , but again not as 
a correlation but in the sense that the line is seen in emission 
at low M em and absorption for high M env . Such behaviour is ex- 
pected if the line optical depth, and thus the H2O column density, 
increases with envelope mass. An optical depth of -unity for 
the Hi s O line is confirme d by detections of HJ 7 Q line emissio n 
in several of our sources (Herpin et al.l 120121 IChoi et al1l2013l) . 
Apparently the average H2O abundance is rather constant across 
all envelopes in our source sample. The studies mentioned in §Q] 
indicate outer envelope abundances of ~10 -9 — 10 8 . 

The increase of the 987 GHz line flux with Lboi is probably 
due to a higher temperature of the inner envelopes of sources 
with a higher luminosity. This idea is confirmed by the trend be- 
tween 987 GHz line flux with far-infrared continuum flux den- 
sity (Fig.fTTTi. Such a trend is seen both at 60 and at 100 //m, and 
like the trend with Lboi, is more a threshold effect than a corre- 
lation. The large beam size of IRAS may influence these trends 
quantitatively, but probably has no strong qualitative effect be- 
cause our sources were selected to be relatively isolated (®. 
Far-infrared pumping could influence the results for individual 
H2O lines, but is unlikely to dominate here, because the effect is 
seen in several lines. 

Further evidence that the H2O 987 GHz line flux is due to en- 
velope heating is presented in Fig. [12] which is a plot of the flux 
in this line versus the mass-weighted temperature from the enve- 
lope models in 34.11 While a 1 : 1 correlation cannot be claimed, 
the points show a clear trend toward higher line fluxes at higher 
average temperatures. 

Figures IPX21 and |D~3"1 show a similar analysis for the narrow 
and broad outflow components in the H2O line profiles of our 
sources. The most significant trend is that the H^O line flux 
of the broad outflow component decreases with increasing Lboi 
(r=-0.683, P-3%), but only in the 11 cases where this com- 
ponent is detected in H^O. Second, the 987 GHz line flux of 
the broad outflow component appears to increase with increas- 
ing M env , but with r=0.174 and P=14%, this trend is less sig- 
nificant than the first one. The line fluxes of the narrow outflow 
component do not show any clear trends with Lboi or M env . 

The lack of evolutionary trends in the H?Q line fluxes of 
our so urces is in contrast with the study of lLopez-Sepulcre et alj 
d2011l) , who observed the SiO 7=2-1 and 3-2 lines toward 57 
high-mass star-forming regions, and found that the SiO line lu- 
minosity drops with evolutionary stage as measured by increas- 
ing values of L\, i/M em . This drop suggests a decline of the jet 
activity with the evolutionary state of the source, a s also ob- 
served for low-mass sources (Bonte mps et al.l [1996). although 
excitation effects also may play a role. 

One reason for this difference may be that even the H2O 
line signals from the outflows may be somewhat optically thick, 
so that saturation masks any evolutionary trends. Therefore, we 
consider next trends of the line width with Lboi and M env . The 
line width is less sensitive to the optical depth than the line in- 
tensity, and has the additional advantage that the appearance of 
the line in emission or absorption does not influence its measure- 
ment. Section 13.41 already shows the trend of the envelope A V 
with M env , which with r=0.407 and P=4% for N=22 is likely 
real. 



1(T 



10 



10 



Clean— up phase 






A 

m 
■ 


• 


• 

A 


• 






< 




• 


• 




* 


■ 




D □ 


■ n 






- 




- 


D .. ' 






Accretion phase 








□ 


1' 






, 



10 



10 



10' 



(M ) 



Fig. 13. Bolometric luminosity versus envelope mass for our 
sources, where the dotted line divides the two phases of high- 
mass protostellar evolution suggested by Molinari et al. (2008) 
and the dashed arrow is a typical evolutionary track. The typi- 
cal uncertainty of a factor of 2 in both luminosity and mass is 
indicated at the bottom. Symbols are as in Fig.fTTI 



Figures lD~4l and 1531 show scatter plots of the widths of the 
narrow and broad outflow components versus Lboi and M env . The 
only significant trend that appears is that AV of the narrow out- 
flow component increases with M env . A relation of this com- 
ponent wit h the envelope has bee n found before for low-mass 
protostars (Kristens en et alj |2012). and suggests that the narrow 
outflow component is due to outflow-envelope interactions. 

4.3. Evolutionary trends 

We now discuss several proposed evolutionary indicators for 
high-mass star formation and test them on our source sample. 
We will use the ratio of Lboi over M env as pri mary evolution- 
ary indica tor, because of previou s success (Van de r Tak et al.l 
2000; iMotte et al.l 120071: lLopez-Sepulcre et al. 2010), because 
it is a continuous variable, and because it has a physical basis 
(declining acc retion activit y and envelope dispersal). In partic- 
ular, Molinari et al.l (120081) have proposed an evolutionary se- 
quence for high-mass protostars in terms of two parameters: 
the envelope mass and the bolometric luminosity. They suggest 
a two-step evolutionary sequence: protostars first accrete mass 
from their envelopes, and later disperse their envelopes by winds 
and/or radiation. In the diagram in Fig. Qj] their calculated evo- 
lutionary tracks start at the bottom, rise upward until they reach 
the dotted line, and then proceed horizontally to the left. The po- 
sitions of our sources in this diagram suggests that the subsam- 
ple of mid-IR-quiet protostellar objects is the youngest, and that 
the other subsamples contain more evolved sources. The same 
effect is seen in Figs. ID.ll4D.3l where the mid-IR-quiet proto- 
stellar objects occupy the left (low-L DO i) sides of the diagrams. 
On the other hand, the presence of a hot molecular core and/or 
an ultracompact HII region is unrelated to mid-IR brightness and 
Lbo\/M em ratio, at least for our sample. The occurrence of these 
phenomena may depend more on absolute luminosity than on 
evolutionary stage. 
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To test the use of H2O line profiles as evolutionary indica- 
tors, we use the appearance of the envelope component in the 
H2O ground-state line profiles, as introduced in 33.41 In Type 1 
sources, all lines appear in absorption; Type 2 sources show 
emission in the 987 GHz line and absorption in the others; in 
Type 3 sources, both the 987 GHz and the H' s O lines appear in 
emission; and Type 4 sources show all lines in emission. While 
the Type 1 sources all seem to be in early evolutionary stages 
as probed by a low Lboi/M env ratio, the other source types are all 
mixed. Thus, an H2O line profile with pure absorption appears 
to be a sign of an early evolutionary stage, whereas emission is 
a sign of a more evolved stage. 

The appearance of the broad and narrow outflow components 
in the H2O ground-state line profiles differs between our sources, 
and for each component, four types may be defined as for the 
envelopes ( {33.51 l. While these emission/absorption types may be 
useful to trace protostellar evolution, our data do not allow to 
test this hypothesis, since most of our sources belong to just one 
of these types, with the broad outflow appearing in absorption 
and the narrow outflow in emission. This issue may be revisited 
after HIFI data become available for a larger sample of sources 
outside the WISH program. 

Besides the shape of the line profile, the fluxes and widths 
of the H2O line may also be evolutionary indicators. The bottom 
panels of Figs. ID. ll4D.3l show scatter plots of our observed line 
fluxes in the envelope and outflow components versus the ratio 
M env /Lboi which as discussed above is sometimes used as evolu- 
tionary indicator for high-mass protostars. No trends appe ar, and 
the same is true for the plots of the line width versus M env /Lboi 
(Figs.l5llLT4landlD31). 

5. Summary and conclusions 

This paper presents spectroscopic observations of low-/ H2O 
transitions observed toward 19 regions of high-mass star for- 
mation. The line profiles show contributions from protostellar 
envelopes, bipolar outflows, and foreground clouds. These con- 
tributions, both in emission and in absorption, are disentangled 
by fitting multiple Gaussians to the line profiles. 

Our observations show that the ground state lines of H2O 
are powerful tracers of the various physical and kinematic com- 
ponents of the high-mass protostellar environment, as well as 
of foreground clouds. The envelope components of the ground 
state lines of H2O at 1113 and 1669 GHz do not show signif- 
icant trends with bolometric luminosity, envelope mass, or far- 
infrared flux density or colour, indicating that these lines trace 
the outer parts of the envelopes where the effects of star forma- 
tion are hardly noticeable. On the other hand, broad wings due to 
outflows are clearly seen in these lines, which are a clear effect 
of star formation. The competing effects of line emission and 
absorption in the ground-state lines lead to low integrated line 
fluxes, which may explain the difficulty in detecting these lines 
in distant galaxies. 

Unlike the ground-state lines, the flux in the excited-state 
987 GHz line increases systematically with luminosity and with 
the far-infrared flux density, both at 60 and at 100 pm. This line 
therefore appears to be a good tracer of the average temperature 
of the source. Indeed, this line and other excited-state lines have 
proven to be more easily detec ted at high redshift (lOmont et alJ 
l20llllVan der Werf et al.ll20lTh . 

Furthermore, the H^O ground-state line turns from emission 
into absorption with increasing envelope mass. This trend indi- 
cates that the H2O abundances in protostellar envelopes do not 
change much with time. 



The appearance of the H2O line profiles in our sample 
changes with evolutionary state, in the sense that the youngest 
sources (as probed by mid-infrared brightness and L\, i/M env ra- 
tio) show pure absorption profiles for all lines, while for more 
evolved sources, one or more of our studied lines appears in 
emission. This effect is visible for the envelope component but 
not for the outflows. In addition, the presence of infall or expan- 
sion, as probed by (inverse) P Cygni profiles, does not appear 
related to evolutionary stage. 

In the future, our team will make detailed estimates of the 
H2O abundances in these sources as a function of radius, based 
on HIFI observations of *20 lines of H 2 0, H' g O and H' 7 
The Herschel/PACS instrument will be used to study higher- 
excitation H2O lines and to improve estimates of the bolometric 
temperatures of our sources, which may act as evolutionary in- 
dicator. Maps of selected lines with HIFI and PACS will clarify 
the spatial distribution of H2O, both on the scale of protostellar 
envelopes and outflows as on larger (protostellar cluster) scales. 
On smaller scales, ALMA will provide images of high-excitation 
H2O and H^O lines and thus zoom into the close vicinity of the 
central protostars. 
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Table A.l. List of observations. 



Source 


H 2 & Hf 


H 2 


H 2 




Li - Ooo 


2(K - Li 


2l2 _ loi 


IRAS 05358 


206126 


204510 


203954 


IRAS 16272 


214419 


203168 


192584 


NGC 6334I(N) 


206385 


204520 


214454 


W43-MM1 


191670 


191616 


192575 


DR21(OH) 


194794 


195026 


192569 


W3 IRS5 


191658 


191612 


192566 


IRAS 18089 


229883 


215911 


216696 


W33A 


191638 


191636 


192576 


IRAS 18151 


229880 


218211 


218909 


AFGL 2591 


197973 


195019 


192570 


G327 


214425 


203170 


192585 


NGC 63341 


206383 


204519 


214455 


G29.96 


191668 


191617 


192574 


G31.41 


191671 


191615 


192573 


G5.89 


229888 


218120 


216694 


G 10.47 


229884 


215914 


216695 


G34.26 


191674 


194995 


192572 


W51e 


194801 


195014 


192571 


NGC 7538 


191663 


201599 


200758 



Notes. The preceding 1342 has been omitted from all obsIDs. 



Appendix A: Observation log and line flux tables 



14 



F.F.S. van der Tak et al.: Low-./ H2O line profiles toward high-mass star-forming regions 
Table A.2. Central velocities, FWHM widths, and integrated line fluxes measured for the protostellar envelopes. 



Source 


^LSR 


AV 


1113 GHz 


H' 8 


1669 GHz 


987 GHz 


IRAS 05358 


-17.9(6) 


3.8(14) 


-13.4 


-0.49 


>17.4 


-1.39 


IRAS 16272 


-45.9(7) 


4.7(20) 


-7.51 


-0.26 


>13.5 


+6.83 


NGC 6334I(N) 


-2.4(5) 


3.9(4) 


-2.96 


-0.70 


>12.3 


-11.8 


W43 MM1 


+ 100.5(7) 


6.0(15) 


-12.9 


-2.07 


-10.1 


-2.17 


DR21(OH) 


-3.3(4) 


3.4(13) 


>12.0 


-0.68 


>9.9 


-0.69 


W3 IRS5 


-37.7(12) 


3.2(6) 


-3.06 


+ 1.84 


-2.00 


+21.9 


IRAS 18089 


+33.7(11) 


3.1(11) 


-3.47 


-0.45 


-8.44 


+0.90 


W33A 


+35.1(12) 


5.6(13) 


-4.34 


-0.61 


>21.0 


<0.60 


CI 


+39.3(6) 


4.8(14) 


-2.85 


+ 1.09 


-5.66 


+9.46 


IRAS 18151 


+33.1(2) 


1.8(10) 


-0.52 


<0.07 


>5.63 


+ 1.35 


AFGL 2591 


-3.6(14) 


4.0(10) 


+24.0 


+0.92 


+8.67 


+ 11.8 


G327-0.6 


-42.5(5) 


3.3(9) 


-1.96 


-0.62 


>15.3 


-1.43 


NGC 63341 


-5.7(16) 


5.7(7) 


-5.31 


-0.79 


-5.46 


+ 15.3 


G29.96 


+98.2(6) 


4.3(13) 


-0.97 


+ 1.10 


>13.8 


+ 19.2 


G31.41 


+99.1(4) 


5.2(14) 


>16.7 


-1.68 


>23.2 


-1.67 


G5.89 


+ 12.0(4) 


7.2(20) 


-5.71 


+4.85 


>22.8 


+70.0 


G 10.47 


+70.0(10) 


3.9(4) 


-1.88 


-0.68 


-5.82 


+4.17 


a 


+65.5(9) 


4.3(13) 


-1.63 


+0.50 


-6.43 


+ 18.3 


G34.26 


+61.5(8) 


4.9(11) 


-22.1 


-0.63 


-8.40 


-1.35 


CI 


+57.1(4) 


5.0(6) 


-5.81 


+2.99 


-7.85 


+27.4 


W51e 


+55.7(17) 


6.3(12) 


-12.5 


-0.81 


-10.2 


+43.8 


NGC 7538 


-57.0(3) 


3.6(10) 


-4.10 


+0.34 


-8.02 


+ 18.7 



Notes. Positive fluxes denote integrated intensities ( fr m bdV in Kkm s _1 ) of emission features, where upper limits denote non-detections. Negative 

fluxes denote integrated optical depths ( frdK in km s _1 ) of absorption features, where lower limits denote saturated absorptions. Central velocities 

and FWHM line widths are in units of kms~', and numbers in parentheses are l<x errors in units of the last decimal. Line fluxes have uncertainties 

of*10%. 

'• a) Source with two envelope components. 
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Table A.3. As TablefO for the outflows. 



Source 


Vlsr 


AV 


1113 GHz 


HfO 


1669 GHz 


987 GHz 


IRAS 05358 


-13.13(54) 


16.9(46) 


+ 14.2 


+0.42 


+ 19.0 


+23.3 


IRAS 16272 


-34.8(3) 


6.72(3) 


-1.73 


n/d 


-6.75 


n/m 




-48.1(28) 


20.4(49) 


+6.16 


n/d 


+3.79 


+ 11.4 


NGC6334I(N) 


-4.3(19) 


9.2(66) 


-3.58 


-2.06 


>14.2 


n/m 




+ 1.0(30) 


17.1(39) 


+21.6 


+4.33 


n/m 


+41.0 


W43 MM1 


+95.85(35) 


6.65(81) 


-8.25 


n/d 


-11.3 


n/m 




+ 105.4(35) 


10.0(48) 


+5.62 


n/d 


+3.41 


+29.9 


DR21(OH) 


-1.15(60) 


2.17(7) 


-2.61 


-0.24 


-2.50 


n/m 




-3.9(46) 


13.1(64) 


+41.0 


+ 1.21 


+21.9 


+80.8 


W3 IRS5 


-41.88(90) 


3.0(14) 


-0.17 


+ 1.74 


^1.57 


+2.85 




-36.6(11) 


21.7(73) 


+26.0 


n/m 


+65.5 


+66.3 


IRAS 18089 


+32.5(27) 


4.6(9) 


-10.8 


+0.59 


>17.0 


n/m 




+37.3(27) 


14.2(28) 


+5.41 


n/m 


+5.67 


+10.5 


W33A 


+29.8(15) 


5.8(26) 


-16.80 


-0.28 


-6.34 


n/m 




+42.6(45) 


20.6(91) 


+ 11.1 


+ 1.09 


+ 10.7 


+ 15.2 


IRAS 18151 


+34.3(4) 


15.6(51) 


+3.98 


n/d 


+4.21 


+4.20 


AFGL 2591 


-10.0(17) 


15.0(16) 


-20.6 


n/d 


-38.3 


+6.30 


G327-0.6 


-48.3(15) 


6.5(26) 


-6.91 


-0.71 


> 18.44 


+3.75 




-37.5(50) 


15.3(47) 


+2.29 


+2.51 


n/m 


+32.6 


NGC 63341 


-10.2(49) 


6.3(23) 


-2.21 


-1.19 


n/m 


-1.22 




-11.0(56) 


14.7(54) 


-1.57 


-0.29 


-21.8 


+55.2 


G29.96 


+ 100.7(56) 


3.9(15) 


-16.96 


n/d 


-2.98 


n/m 




+99.13(39) 


18.4(24) 


+ 16.1 


n/d 


+24.6 


+20.2 


G31.41 


+95.0(4) 


4.0(19) 


-5.36 


n/d 


-2.34 


+6.91 




+99.65(3) 


34.0(3) 


n/m 


n/d 


n/m 


+ 12.8 


G5.89 


+44.1(30) 


14.8(4) 


+7.55 


n/m 


+ 18.46 


n/m 




+27.7(32) 


19.3(65) 


+20.5 


+3.86 


+35.71 


+ 119.4 




-9.4(28) 


30(13) 


-3.78 


-1.64 


+34.14 


+65.7 




-26.9(5) 


13.9(45) 


-0.70 


-0.62 


n/m 


n/m 




-49.5(85) 


16.2(72) 


-1.14 


-0.30 


-1.53 


+4.83 


G 10.47 


+62.0(5) 


4.1(5) 


+2.70 


n/m 


n/m 


+4.45 




+80.2(51) 


8.3(6) 


+4.71 


n/m 


-1.36 


+5.47 


G34.26 


+53.7(77) 


4.2(13) 


-9.60 


+ 1.20 


-9.14 


n/m 




+57.0(85) 


18.1(81) 


-1.66 


-0.96 


-10.7 


+25.8 


W51e 


+67.7(13) 


6.9(28) 


-14.3 


-0.25 


-28.2 


-1.44 




+61.7(13) 


14(11) 


-17.6 


-0.91 


^t.70 


+62.5 


NGC 7538 IRS 1 


-57.7(51) 


2.25(43) 


-2.36 


-0.13 


-0.80 


n/m 




-56.8(58) 


13.4(24) 


+ 10.3 


+ 1.67 


+2.66 


+ 17.4 



Notes. The symbol n/d means that the outflow is not at all detected in this line; the symbol n/m means that only this component has no match in 
this line. 
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Table A.4. As Table lA.2l for the foreground clouds. 



Source 



V,: 



AV 



1113 GHz 1669 GHz N(H 2 0) o/p 



IRAS 16272 


-41.38(8) 


1.67(3) 


-2.16 


>6.33 


>35 


>5.9 




-38.80(23) 


1.94(5) 


-1.36 


-2.50 


15 


3.7 


NGC6334I(N) 


+6.86(13) 


1.49(12) 


-1.52 


>5.12 


>27 


>6.6 


W43 MM1 


+62.1(6) 


1.04(8) 


-0.094 


-0.11 


0.73 


2.4 




+67.21(4) 


1.81(8) 


-1.13 


-1.80 


11 


3.2 




+70.96(8) 


1.6(2) 


-0.73 


-1.15 


7.1 


3.2 




+77.(7) 


4.1(6) 


-1.63 


-1.11 


9.0 


1.4 




+79.7(1) 


1.93(5) 


-1.48 


-3.67 


21 


5.0 




+82.4(2) 


2.49(5) 


-1.63 


-2.58 


16 


3.2 




+87.83(12) 


2.1(4) 


-0.20 


-0.49 


2.75 


4.9 


DR21(OH) 


+ 14.9(1) 


1.5(2) 


-0.35 


-0.90 


5.0 


5.1 




+12.3(1) 


1.8(2) 


-0.47 


-2.03 


11 


8.6 




+ 10.0(1) 


1.7(7) 


-0.51 


-2.53 


13 


9.9 




+8.4(4) 


1.9(2) 


-3.78 


-0.80 


13 


0.4 




+6.9(1) 


1.22(1) 


-0.51 


-0.38 


3.0 


1.5 


W3 IRS5 


-20.8(1) 


1.04(8) 


-0.61 


>3.56 


>18 


>12 


IRAS 18089 


+17.6(1) 


3.4(6) 


-2.81 


-7.24 


40.3 


5.2 




+ 19.7(3) 


2.8(11) 


-0.91 


-2.34 


13.0 


5.1 




+21.1(1) 


2.2(4) 


-0.95 


-1.42 


8.8 


3.0 




+23.9(2) 


1.9(1) 


-0.82 


-2.68 


14.4 


6.5 




+25.7(2) 


2.67(7) 


-1.46 


^.00 


22.0 


5.5 


W33A 


+24.01(7) 


1.5(1) 


-0.66 


-1.26 


7.4 


3.8 


IRAS 18151 


+27.7(2) 


0.7(2) 


-0.36 


-1.06 


5.8 


5.9 


AFGL 2591 


-0.67(13) 

+ 13° 
-38.02(4) 


3.0(6) 


>9.54 


>7.63 


>58 


-1.6 


G327-0.6 rf 


1.18(33) 


-0.06 


-0.19 


1.1 


6.0 


NGC 63341 


+0.48(11) 


2.58(52) 


-1.45 


-3.42 


13.8 


2.7 




+6.63(17) 


1.61(16) 


-2.11 


>5.63 


>31.2 


>5.3 




+8.72(39) 


1.27(52) 


-0.20 


-0.17 


1.3 


1.7 


G29.96 


+ 103.5(1) 


1.7(1) 


-1.09 


-2.91 


16 


5.3 




+91.6(2) 


3.4(5) 


-0.67 


-0.59 


4.3 


1.8 


G31.41 


+83.0(3) 


3.4(6) 


-0.48 


-1.27 


7.0 


5.3 


G5.89 


+ 19.85(1) 


4.4(6) 


-1.90 


-25.17* 


121.7* 


26.5* 


G 10.47 


+42.3(1) 


1.1(3) 


-0.16 


-0.24 


1.5 


3.0 




+47.06(4) 


1.35(7) 


-0.09 


-0.19 


1.1 


4.5 




+84.9(1) 


2.6(10) 


-0.19 


-1.67* 


8.2* 


17.7* 




+91.21(6) 


1.4(4) 


-0.08 


-0.32 


1.7 


8.3 


G34.26 


+11.4(1) 


2.6(1) 


-6.43 


-14.0 


80 


4.4 




+ 14.26(3) 


1.18(1) 


-0.35 


-0.55 


3.4 


3.1 




+27.09(1) 


1.6(3) 


-1.52 


-5.48 


29 


7.2 


W51e 


+48.66(1) 


2.27(1) 


-0.42 


-1.14 


6.3 


5.4 




+45.57(9) 


1.2(1) 


-2.51 


>4.42 


>26 


>3.5 




+ 12.7 C 














+6.48(5) 


2.4(2) 


-2.63 


-2.99 


20 


2.3 




+4.87(1) 


1.1(2) 


-0.27 


-0.84 


4.5 


6.2 



Notes. The last two columns list the derived HtO column density (in 1 12 cm -2 ) and orth o/para ratio. 

(fl) Blended with H 2 1661 GHz line emission from other sideband; see lChoi et al.l ( 120130 

'*' Highly uncertain because 1669 GHz line is saturated. 

(c) Only detected at 1 1 13 GHz. 

{d} WBS spectra of this source show spatially extended H2O foreground absorbers between -20 and -3 km s (S. Leurini, priv. comm.). 
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Appendix B: Observed line profiles 




Fig.B.l. Spectra of the H 2 In - o (bottom), H' 8 In - o 
(second from bottom), H2O 2i2-loi (third from bottom) and 
H 2 2 2-lii (top) lines toward IRAS 05358. The vertical scale 
is T m b in K, and the dashed line denotes the systemic velocity 
determined from ground-based observations. The spectra have 
been shifted vertically for clarity. 
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Fig.B.2. As FigurelO toward IRAS 16272. 
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Fig.B.3. As FigureEU toward NGC 6334I(N). 
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Fig.B.4. As FigureEU toward W43 MM1. 
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Fig.B.5. As FigureED toward DR21(OH). 
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Fig.B.6. As FigureEU toward W3 IRS5. 
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Fig.B.7. As FigurelO toward IRAS 18089. 
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Fig. B.8. As Figure lETTl toward W33A. The emission/absorption 
feature near km s" ' in the 1 669 GHz spectrum is the H2O 1 66 1 
GHz line from the other sideband. 
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Fig.B.9. As FigurelO toward IRAS 18151. 
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Fig.B.ll. As Figure EU toward G 327. 
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Fig. B.10. As FigureEm toward AFGL 259 1 . The emission fea- . . 

ture near +18 km s" 1 in the 1669 GHz spectrum is the H 2 1661 Fl §- BJ2 - As FigureED toward NGC 63341. 
GHz line from the other sideband. 
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Fig.B.13. As Figure|RT] toward G29.96. 
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Fig.B.14. As Figure |BJ] toward G31.41 




Fig. B. 15. As Figure |BT1 toward G5.89. The emission feature 
near +90 kms" 1 in the H^O spectrum is the 13 CO 10-9 line. 
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Fig.B.16. As FigureET] toward G10.47. 
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Fig.B.17. As FigurelO toward G34.26. 
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Fig. B.19. As Figure lO toward NGC 7538 IRS 1 . 



Fig.B.18. As FigurelO toward W51e. 
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Appendix C: Continuum models 
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Fig. C.l. Continuum model for IRAS 05358. Top left: Spectral energy distribution. Top right: Submillimeter image with 3<x contour 
marked in white and useable map sectors marked in black. Bottom left: Temperature and density structure as a function of radius. 
Bottom right: Submillimeter emission profiles. The numbers in the top left panel are the modeled envelope mass, the modeled 
luminosity, the adopted envelope size and the adopted distance. 
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Fig. C.2. As FigureO for IRAS 16272. 
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Fig. C.3. As FigureO for NGC 6334I(N). 
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Fig. C.4. As FigureO for W43-MM1. 
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Fig. C.5. As Figure O for DR21(OH). 
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Fig. C.6. As Figure O for W3 IRS5. 



27 



F.F.S. van der Tak et al.: Low-./ H2O line profiles toward high-mass star-forming regions 



10 r 172 M„ 
1.4e4 L„ 
3.9e4 AU 
„ 2.3 kpc 

10 3 r 



IRAS 18089-1732 




10 10 

Wavelength [j(im] 




56.0 18:11:52.0 48.0 

Right ascension 



io'L 
10 s - 



p = 1.1 



10 10 10 

Radius [AU] 




ID 1 






N 


















CO 






a 












& 0.1 




%^ ] 




850 ^m 


W 1 ^ : 


0.01 





10 
Offset [arcsec] 



Fig. C.7. As FigureO for IRAS 18089. 
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Fig. C.8. As Figure|CT] for W33A. 
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Fig.C.9. As Figure O for IRAS 18151. 
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Fig. CIO. As FigureO for AFGL 2591. 
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Fig. C.ll. As FigurelCU for G327-0.6 
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Fig. C.12. As FigureO for NGC 63341. 
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Fig. C.13. As FigureO for G29.96. 
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Fig. C.14. As FigureO for G3 1 .41 . 
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Fig. C.15. As Figurelcn for G5.89 
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Fig. C.16. As FigureO for G10.47. 
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Fig. C.17. As FigureO for G34.26. 





10 10 

Wavelength [/zm] 



44.0 19:23:40.0 

Right ascension 




10 10 10 

Radius [AU] 



Fig. C.18. As FigureO for W51e. 
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Fig. C.19. As FigureO for NGC 7538 IRS1. 
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Fig. D.l. Observed line fluxes (positive: emission) or ab- 
sorbances (negative: absorption) versus bolometric luminosity L 
of the source (top), versus envelope mass M (middle), and versus 
the ratio L/M (bottom). Open squares: mid-IR-quiet HMPOs; 
filled squares: mid-IR-bright HMPOs; triangles: hot molecular 
cores; circles: ultracompact HII regions. Uncertainties in line 
fluxes are a:10%, while the typical uncertainty of masses, lumi- 
nosities and L/M ratios is a factor of 2. 
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Fig. D.2. As Figure lETTl for the narrow outflow component. 



Fig. D.3. As Figure lDTI for the broad outflow component. 
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Fig. D.4. Line width of the narrow outflow component versus Fig. D.5. As Figure lLVfl for the broad outflow component, 
bolometric luminosity L of the source (top), versus envelope 
mass M (middle), and versus the ratio L/M (bottom). The er- 
ror bar denotes the typical uncertainty in mass, luminosity, and 
L/M ratio; the uncertainty in line width is smaller than the sym- 
bol size. 
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